ON THE BEHAVIOUR OF A RUMOUR PROCESS 
WITH RANDOM STIFLING 



ELCIO LEBENSZTAYN, FABIO P. MACHADO, AND PABLO M. RODRIGUEZ 

Abstract. We propose a realistic generalization of the Maki-Thompson rumour model 
by assuming that each spreader ceases to propagate the rumour right after being involved 
in a random number of stifling experiences. We consider the process with a general initial 
configuration and establish the asymptotic behaviour (and its fluctuation) of the ultimate 
proportion of ignorants as the population size grows to oo. Our approach leads to explicit 
formulas so that the limiting proportion of ignorants and its variance can be computed. 



1. Introduction 

In the past decades, there has been great interest in understanding and modelhng differ- 
ent processes for information diffusion in a population. Most of the time, the mathematical 
theory of epidemics is adapted for this purpose, even though there are differences between 
the process of spreading information and the process of spreading a virus or a disease. 
In the standard versions of the models, the most noticeable differences are between the 
way spreaders cease to spread an item of information and the way infected individuals are 
removed from epidemic processes. Still, some slightly modified models fit both processes 
(see, for example, Dunstan [7], where the general stochastic epidemic model is considered 
as a model for the diffusion of rumours). 

Kurtz et al. [l^ recently introduced a model using a complete graph in which, as soon 
as an individual is infected, an anti- virus is given to that individual in such a way that the 
next time a virus tries to infect it, the virus is ineffective. Besides, a virus can survive up 
to L individuals empowered with anti-virus. Individuals are represented by the vertices 
of the complete graph, while the virus is represented by a moving agent that replicates 
every time it hits a healthy individual. The authors prove a Weak Law of Large Numbers 
and a Central Limit Theorem for the proportion of infected individuals after the process 
is completed. 

There are two classical models for the spreading of a rumour in a population, which were 
formulated by Daley and Kendall Q and Maki and Thompson [l^. In the model proposed 
by Maki and Thompson [l3|, a closed homogeneously mixing population experiences a 
rumour process. Three classes of individuals are considered: ignorants, spreaders and 
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stiflers. The rumour is propagated through the population by directed contact between 
spreaders and other individuals, which are governed by the following set of rules. When a 
spreader interacts with an ignorant, the ignorant becomes a spreader; whenever a spreader 
contacts a stifier, the spreader turns into a stifier and when a spreader meets another 
spreader, the initiating spreader becomes a stifier. In the last two cases, it is said that the 
spreader was involved in a stifling experience. Observe that the process eventually ends 
(when no more spreaders are left in the population). 

We show how the techniques used by Kurtz et al. [lo| in the context of epidemic models 
can be useful in studying a general rumour process. In particular, we propose a generaliza- 
tion of the Maki- Thompson model. In our model, each spreader decides to stop propagating 
the rumour right after being involved in a random number of stifling experiences. 

To define the process, consider a closed homogeneously mixing population of size + 1. 
Let i? be a nonnegative integer valued random variable with distribution given hj P{R = 
i) = Tj for i = 0, 1, . . . , and let fi = E[R] > and = Var[_R]. Assign independently 
to each initially ignorant individual a random variable with the same distribution as R. 
Once an ignorant hears the rumour, the value of R assigned to him determines the number 
of stifling experiences the new spreader will have until he stops propagating the rumour. 
If this random variable equals zero, then the ignorant joins the stiflers immediately after 
hearing the rumour. 

For i = 1,2, ... , we say that a spreader is of type i if this individual has exactly i 
remaining stifling experiences. We denote the number of ignorants, spreaders of type i and 
stiflers at time t by X^^\t), Yi^^\t) and Z^^\t), respectively. Let Y^^\t) = Yi^^\t) 
be the total number of spreaders at time t, so X^^^ (t) + Y^^^ {t) + Z(^) {t) = N + 1 for all t. 
Notice that the infinite-dimensional process 

{V^'^Hm.^o ■■= {(xW(t),FrW,>^f Ht),...)}t>o (1.1) 
is a continuous time Markov chain with increments and corresponding rates given by 

increment rate 
(-1,0,0,...) roXY 

(-l,0,...,V,i,T,...) nXY 

(0,...,0,'l\-l,T,...) {N-X)Y, 
(0,-1,0,0,...) {N-X)Y,. 

We see that the first case indicates the transition of the process in which a spreader 
interacts with an ignorant and the ignorant becomes a stifier immediately (which happens 
with probability tq). The second case indicates the transition in which a spreader interacts 
with an ignorant and the ignorant becomes a spreader of type i (which happens with 
probability r^). The third case represents the situation in which a spreader of type i is 
involved in a stifling experience but remains a spreader (of type i — 1), and finally the last 



i = 1,2,... 
1 = 2,3,... 
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transition indicates the event that a spreader of type 1 is involved in a stifling experience, 
thus becoming a stifler. 

We suppose that the process starts with 

X(^)(0) = (iV+l)xr\ 

f/^)(0) = (iV + 1) y^^ for ^ = 1, 2, . . . and 

zW(0) = (iV+l)4^). 

That is, x\^\yfl\ z\^'^ G [0, 1] are the initial proportions of ignorants, spreaders of type i 
and stiflers of the population, respectively, which are defined in such a way that 

oo 

EC>0 

i=l 

oo 
i=l 

In addition, we assume that the following limits exist: 

Xq = lim xi^^ > and 

A''— >-oo 

.{N) 

L 

and that 



yi,o = lim ylo' for alH = 1, 2, . . . , 



oo 



Wo = ^iyifl < oo. 



i=l 



As already mentioned, the process eventually ends. Let 

r(^)=inf{t:r(^)(t) = 0} 

be the absorption time of the process. Our main purpose is to establish limit theorems 
for the proportion X^^\t^^^) of ignorants at the end of the process. For the classical 
Maki-Thompson model, this problem was first studied rigorously by Sudbury p^], who 
proved, by using martingale arguments, that this proportion converges in probability to 
0.203. This result was later generalized by Watson 1^ using the normal asymptotic 



approximation. Lefevre and Picard [12| derived the exact joint distribution of the final 
number of people who heard the rumour and the total personal time units during which 
the rumour was spread. In Belen and Pearce 0], the authors present an analysis of the 
proportion of the population who never hear the rumour starting from a general initial 
condition. See also Chapter 5 of Daley and Gani [s*] for an excellent account of rumour 
models. 

The approach used to prove our theorems is the theory of density dependent Markov 
chains, presented in Ethier and Kurtz [8]. To the best of our knowledge, this technique is 
used for the first time in the context of rumour models here and in Lebensztayn et al. [ll[ . 
In that paper, the authors study a family of rumour processes which includes the classical 
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Daley-Kendall and Maki-Tliompson models as particular cases. The results presented here 
are of independent interest, as they refer to a generalization of the Maki- Thompson model 
with random stifling and general initial configuration. 

2. Main results 

Definition 2.1. Suppose that /i < oo and consider the function / : (0,xo] — > M given by 

X 

f{x) = Wo + (1 + yu)(xo - x) + log — . 

We define x^o = Xoo{fi,Xo,Wo) as the unique root of / in the interval (0,Xo] satisfying 
fix) > 0. 

Notice that Xqo is the unique root of /, except in the case where Xq > (1 + /i)"^ and 
Wq = 0. See Figure [TJ 



.To > 1/(1 + li) and Wq > Xo < 1/(1 + fi) and tuo > 





/ T 


/ ■ X 





/.r^. .To 










.r„ > 1/(1 + /j) and tui, = Tq < 1/(1 + A*) and iii„ = 







T. 















Figure 1. Behaviour of / - The four possible cases in terms of Xq and Wq. 



Remark 2.2. We can express x^o in terms of the Lambert W function, which is the inverse 
of the function x i— )■ xe^. Indeed, x^o satisfies 

T = T p-(l+M)(^o-3;oo)-i«o 
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which can be written as 

- Xo(l + /i) e-^«(i+^)~"'« = -Xoo(l + /i) e~^°°(^+^'\ (2.1) 

Then, if Wq denotes the principal branch of the Lambert W function (that is, the branch 
that satisfies W{x) > —1), we obtain from f l2.ip that 

Xoo(/i, xo, Wo) = -(1 + fi)-' Woi-xo (1 + fi) e--'o(i+'^)-o), (2.2) 

by noting that — e^^ < — xo(l + fi) e~^o(^+^)~"'o < 0. More details about the Lambert 
function can be found in Corless et al. 

Next, we state the Weak Law of Large Numbers for the proportion of the population 
who have never heard the rumour. 

Theorem 2.3. //O < /i < oo, then 

XW(rW) 
lim — = Xoo in probability. 

N^oo N 

As a consequence of this theorem. 
Corollary 2.4. If jj, = oo, then 

XW(r(^)) 
lim = in probability. 

Proof. Let Ri and R2 be nonnegative integer valued random variables such that Ri < R2 
stochastically, that is, P{Ri > i) < P{R2 > i) for all i > 1. Consider the processes 
{Vl''\t)}t>o and {Vf\t)}t>0 defined as in (11.11) by using the random variables Ri and R2, 
respectively, and with the same initial conditions. Let r{^'' and Tg^'* be the respective 
absorption times. By a coupling argument, these processes can be constructed in such a 
way that 

Xf)(rf))<xr(rr), a.s. (2.3) 

Now suppose that {V^^\t)}t>o is the process defined in (11.11) with the random variable R 
satisfying n = 00. Recall that = P{R = i), i > 0, and for each k > 1 define the random 
variable Rk with distribution given by 



P(i?fc = i) = ri ii i<k and P{Rk = k) = Y^ 

j=k 



By construction, we have that Rk ^ R stochastically for all k. Taking this and (12.31) into 
account, we conclude that 

XW(rW)<xf)(rf)), a.s. (2.4) 
for all k. Then (12.41) and Theorem 12.31 imply that 

XW(rW) 

0<limsup < Xoo{^k,xo,wo) a.s.. 
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where Hk = E[Rk]. Since xq > and limfc_yoo/ifc = oo, we have that xq > (1 + fJ'k)^^ for 
large enough k and in this case a:oo(/^fc, a;o, wo) (given by (12.21) with /i = /ifc) goes to as 
— !■ oo. This completes the proof of CoroUarv 12.41 □ 

We now present the Central Limit Theorem for the ultimate proportion of ignorants in 
the population. 

Theorem 2.5. Suppose that u"^ < oo. Assume also that wq > or that wq = and 
Xq > (1 + Then, 



^/XW(rW) 



N 



N{0,a'^) as N ^oo, 



where =^ denotes convergence in distribution, and N{0, a^) is the Gaussian distribution 
with mean zero and variance given by 



Remark 2.6. Observe that our results refer to a general initial condition, similar to that 
considered in the deterministic analysis presented in Belen and Pearce The process 
starting with one spreader and ignorants corresponds to xq = 1 and wq = 0, in which case 
the limiting fraction of ignorants and the variance of the asymptotic normal distribution 
in the CLT reduce respectively to 

Xoo = Xoo(/i, 1, 0) = -ifi + 1)-' Wo{-{fi + 1) e-(^+i)) and 

2 _ Xoo(l - Xoo)(l - (1 + - Z/^) Xqo) 
~ (l-(l + /i)Xoo)2 

The behaviour of function of n is shown in Figure [21 



,(1 - (Xq ^ + Wo + {xq - Xoo){l + yU - l^^)) Xoo) 
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Figure 2. Graph of Xoo(/i, 1, 0). 
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Here are some important cases: 

(a) For R = K, {k, > 1 an integer), we have the K-fold stifling Maki- Thompson model (so 
called by Daley and Gani jsj, in the context of the Daley-Kendall model), for which 

2 -^ooll -^oo) 
— 7 ^ ; 

1 - + l)a;oo 

where 

Xoo = a:oo(«:, 1, 0) = -{k + 1)-^ Woi^iK + 1) e-^^+'^). (2.6) 

Table [U exhibits the values of Xoo and in this case for k = 1, . . . , 8. The original Maki- 
Thompson model is obtained by considering R = 1, xq = 1 and wq = 0, consequently 
our theorems generalize classical results proved by Sudbury [15] and Watson [l6|. For the 
2-fold stifling Maki- Thompson model, the asymptotic value of x^o ~ 0.0595 was originally 
obtained by Carnal jsj. Formula fl2.6p is presented in Appendix D of Belen's doctorate 
thesis [l!]. 

(b) Let R ~ Geometric(p), that is, ro = and ri = p (1— p)*~^, i = 1,2, . . . In this model, 
an ignorant always becomes a spreader upon hearing the rumour and each time a spreader 
meets another spreader or a stifler, he decides with probability p to become a stifler, 
independently for each spreader and each meeting. Thus, given that a spreader has not 
yet stopped propagating the rumour, the conditional distribution of the additional number 
of stifling experiences he will have does not depend on how many stifling experiences he 
already had. This means that every time a spreader chooses whether or not to become 
a stifler, he does not have a "memory" of how many unsuccessful telling meetings he has 
been involved in. Table [2] shows the values of Xoo and cr^ for a;o = 1, u^o = and some 
arbitrarily chosen values of p. 

(c) Consider R ~ Poisson(A), in which case an ignorant individual has the choice (with 
a positive probability equal to e~'^) of becoming a stifler as soon as he learns the rumour. 
Moreover, in his successive decisions about stifling, a spreader does have some "memory" 
of the number of his previous stifling experiences. Table [3] presents the values of Xoo and cr^ 
for Xq = 1, Wq = and some values of A. 



3. Proofs 

Here are the main ideas in the proofs of Theorems 12.31 and 12.51 First, by means of a 
suitable time change of the process, we define a new process {V^'^\t)}t>Q with the same 
transitions as {V^^^ {t)}t>o, so that they end at the same point of the state space. Next, we 
work with a reduced Markov chain obtained from {V^^\t)}t>Q in order to apply the theory 
of density dependent Markov chains presented in Ethier and Kurtz ^ . As the arguments 
follow a path similar to that presented in Kurtz et al. [lo|, we present only a brief sketch 
of the proofs. 
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K 


1 


2 


3 


4 


5 


6 


7 


8 




0.203 
0.273 


0.0595 
0.0681 


0.0198 
0.0211 


0.00698 
0.00718 


0.00252 
0.00255 


0.000918 
0.000923 


0.000336 
0.000337 


0.000124 
0.000124 



Table 1. /t-fold stifling model, xq = 1 and wq = 0. 



p 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 0.7 


0.8 


0.9 




0.0000167 
0.0000167 


0.00252 
0.00268 


0.0139 
0.0163 


0.0340 
0.0427 


0.0595 
0.0780 


0.0878 0.117 
0.118 0.159 


0.147 
0.199 


0.175 
0.238 



Table 2. i? ~ Geometric (p), xq = 1 and wo = 0. 



A 


0.1 


0.3 


0.5 


0.7 


0.9 


1.1 


1.3 


1.5 


1.7 


1.9 




0.824 
2.908 


0.577 
1.654 


0.417 
1.012 


0.309 
0.654 


0.233 
0.440 


0.178 
0.307 


0.138 
0.219 


0.107 
0.160 


0.0844 
0.119 


0.0668 
0.0895 



Table 3. i? ~ Poisson(A), xq = 1 and wq = 0. 



Time-changed process. Since the distribution of X^^\t^'^^) depends on the process 
{V^^^ {t)}t>Q only through the embedded Markov chain, we consider a time-changed version 
of the process. Let {V^'^\t)}t>o be the infinite-dimensional continuous time Markov chain 

{(x(^)(t),yrw,i^r(t),...)h>o 

with increments and corresponding rates given by 
increment rate 
(-1,0,0,...) roX 

(-i,o,...,V,i,'o\...) nx i = i,2,... 

(0, 0, ... , 'l\ -1, 'o\ . . . ) {N-X)Y, z = 2,3,... 

(0,-1,0,...) {N~X)Y, {¥)-'. 

Furthermore, {V^*-^''(t)}f>o can be defined in such a way that it has the same initial state 
and the same transitions as {V^'^\t)}t>o, so both have the same embedded Markov chain. 
Thus, by defining 

fW = inf{t:F(^)(t)=0}, 
we have that X(^)(r(^)) = XW(fW). 

Dimension reduction and deterministic limit. In order to prove the desired limit 
theorems using Theorem 11.2.1 of Ethier and Kurtz [^], we work with a reduced Markov 
chain. We define 
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and note that the process {{X^^\t), W'^^\t))}t>o is a continuous time Markov chain with 
increments and rates given by 

increment rate 

i^ = {-l,^) TiX 1 = 0,1,... (3.1) 
£„i = (0,-l) N-X. 

Now we define, for t > 0, 

i)W(t) = = iv-i(x(^)(t),iyW(t)), 

and consider 

f3e_j^{x,w) = 1 — X and (3e.{x,w) = riX, i = 0,l,... 
Notice that the rates in (13. ip can be written as 



X W_ 

iV' iV 



so {v^^\t)}t>Q is a density dependent Markov chain with possible transitions in the set 

{U,eo,ii,.r.}. 

Now we use Theorem 11.2.1 of Ethier and Kurtz \&\ to conclude that the process 
{v^^\t)}t>o converges almost surely as — ?■ oo to a deterministic limit. The drift function 
defined in Ethier and Kurtz by F{x, w) = Yli^~i P^ii^i this case given by 

F{x, w) = {—X, (/i + l)a; — 1). 

Hence the limiting deterministic system is governed by the following system of ordinary 
differential equations 

'x'{t) = -x{t), 
w'{t) = {fi + l)x-l 

with initial conditions x(0) = Xq and w{0) = Wq. The solution of this system is given by 
v(t) = {x(t),w(t)), where 

x{t) = Xq e~* and w(t) = f{x{t)) = wq + {1 + fi) {xq — x(t)) — t. 

According to Theorem 11.2.1 of Ethier and Kurtz [Sj, we have that on a suitable prob- 
ability space, 

lim v^^\t) = v{t) a.s. 

iV->-oo 

uniformly on bounded time intervals. In particular, it can be proved that 

lim S:^^\t) = x{t) a.s. (3.2) 

uniformly on M. See Lemma 3.6 in Kurtz et al. [13] for an analogous detailed proof. 
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Proofs of Theorems 12.31 and 12.51 To prove both theorems, we use Theorem 11.4.1 
of Ethier and Kurtz j8|. We adopt their notations, except for the Gaussian process V 
defined on p. 458, that we would rather denote hj U = {Ux, Uw). Here, ip{x,w) = w, and 

Too = inf{t : w{t) < 0} = + (1 + fi){xo - Xoo). 

Moreover, 

Vifiivir^)) ■ F{viT^)) = w'ir^) = (fi + l)x^ - 1< 0. (3.3) 
Proof of Theorem \ 2.3i . We note that wq > and (13. 3p imply that w{too — e) > and 
w{too + e) < for < £ < Too- Then, the almost sure convergence of w^^^ to w uniformly 
on bounded intervals yields that 

lim f(^) = Too a.s. (3.4) 

In the case where wq = and a;o > (1 + /i)^^, this result is also valid because w'{0) > and 
(13. 3p still holds. On the other hand, if i^o = and Xq < {1 + then w{t) < for all 
t > 0, and again the almost sure convergence of w^^^ to w uniformly on bounded intervals 
yields that lim^v-^oo r^^) = = Too almost surely. Therefore, as X(^)(r(^)) = X(^)(f(^)), 
we obtain Theorem 12.31 from (13. 2p and (13.40 . 

Proof of Theorem \2.5[ From Theorem 11.4.1 of Ethier and Kurtz [sj, we have that if > 
or Wo = and Xq > {1 + then y/N {x^^^f'^^'') — Xoo) converges in distribution as 

oo to 

f/x(roo) + , , UUtoo). (3.5) 

(/i + l)Xoo - 1 

The resulting normal distribution has mean zero, so, to complete the proof of Theorem 
12. 5[ we need to calculate the corresponding variance. 

To this end, we have to compute the covariance matrix Cov(t/(roo), U{too)), a task that 
can be accomplished using a mathematical software. The first step is to calculate the 
matrix of partial derivatives of the drift function F and the matrix G. We obtain 

-1 O' 



dF{x, w 
and 

G{x, w) = 



(/i + 1) 



■3j ^Aj'JC 

fix [y^ + /^^ — l)a; + 1^ 
Next, we compute the solution $ of the matrix equation 

^ $(t, s) = dF{x{t),w{t)) $(t, s), $(s, s) = /2, 

which is given by 



V(/^ + l)(l-e-(*-^)) 1^ 
Hence, the covariance matrix of the Gaussian process U at time t is obtained by the formula 

Cov{U{t),U{t)) = [ <^{t,s)G{x{s),w{s))[<!>{t,s)f ds. (3.6) 
Jo 
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As the final step to compute Cov([/(roo), f/(Too)), we have to replace e * and t in the 
formula obtained from (13.61) by Xoo/xq and Too, respectively. The resulting formulas are 

Var([/^(roo)) = ((xq - a;oo)a;oo)/a;o, 

Var(f/^(roo)) = (yU + l)^(a;o - a;oo)a;oo/a;o + ^'^{xq - x^o) 

+ (1 - 2{fi + l)xoo)roo, 

Cov(?7^(roo), f/^(T-oo)) = r^Xoo - (/i + l)(a;o - Xoo)a;oo/xo. 

Using that Too = + (1 + /x)(xo — Xoo), (13.511 and well-known properties of the variance, 
we get formula (12. 5p . 

4. Concluding remarks 

We have proposed a general Maki-Thompson model in which an ignorant individual 
is allowed to have a random number of stifling experiences once he is told the rumour. 
The assigned numbers of stifling experiences are independent and identically distributed 
random variables with mean fi and variance z/^. We prove that the ultimate proportion 
of ignorants converges in probability to an asymptotic value as the population size tends 
to oo. A Central Limit Theorem describing the magnitude of the random fluctuations 
around this limiting value is also derived. The asymptotic value and the variance of the 
Gaussian distribution in the CLT are functions of /i, u"^ and some constants related to the 
initial state of the process. 

We observe that in fact it is possible to obtain another result, concerning the mean 
number m^^^ of transitions that the process makes until ab sorp tion. Using an argument 
analogous to that presented in Theorem 2.5 of Kurtz et al. it can be proved that, if 
z/^ < oo, then 

lim iV~^ m^^^ = Too = u^o + (1 + Ai)(xo - Xoo)- 

A^— >-oo 

As a final remark, we would like to point out the usefulness of the theory of density 
dependent Markov chains as a tool for studying the limiting behaviour of stochastic rumour 
processes. This approach constitutes an alternative to the pgf method and the Laplace 
transform presented in Daley and Gani Gani ^ and Pearce 14 . 
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